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Charm (and Beauty) Production at the Tevatron 

Jonas Rademacker on behalf of the CDF and the D0 Collaboration 
University of Bristol, H H Wills Physics Laboratory, Bristol, BS8 1 TL, UK 

We present recent results on heavy flavour production at Tevatron Run II for typically ~ 1 fb~ 
of analysed pp data at y/s = 1.96 TeV. This includes results on single and correlated open charm 
and bottom cross sections, charm pair production kinematics, J/tp, ^(2S) and XcJ cross sections 
and polarisation measurements in J/f/', ^>(2S), T(1S), and T(2S). 



o 
o 

(N 

> 

O 

oo 

(N 



X 

<D ' 
i ' 

Or 
43. 



> 

in 



o 



X 



INTRODUCTION 



The large bb and cc cross section at s/s = 1.96 TeV 
pp makes the Tevatron a unique place for the study 
of the production and decay of heavy flavour. Al- 
ready in Run I, heavy flavour production measure- 
ments challenged theory, finding heavy flavour pro- 
duction cross sections significantly higher than pre- 
dicted by Next-to-Leading-Order (NLO) QCD calcu- 
lations (see for example [Ij and Q). The discrepancy 
between data and theory was particularly dramatic in 
the quarkonium production, where the "Colour Sin- 
glet Model" leading-order QCD calculation underes- 
timates the measured cross section by more than an 
order of magnitude (see e.g.Q). 



Since the first hadroproduction measurements in 
Run I, a number of theoretical advances have been 
made. Fixed-Order Next to Leading Logarithm 
(FONLL) calculations [J], describe the open charm 
and b production cross sections well. Competing mod- 
els have been put forward that describe the observed 
quarkonium production rates and px spectra well, but 
disagree on their results for quarkonium polarisation. 



There has also been dramatic experimental 
progress. In Run II, which started in 2001, the Teva- 
tron collides p and p at unprecedented luminosity and 
energy, and the D0 and the CDF experiment have 
undergone significant upgrades, many of them opti- 
mising the detector for flavour physics. In this paper, 
we summarise the heavy flavour production measure- 
ments in Run II at D0 and CDF, and compare them 
with theoretical predictions. 

As we will see below, heavy flavour production at 
hadron colliders is a vibrant field which is lead by ex- 
periment rather than theory. A particular challenge 
for theory is the quarkonium polarisation, for which 
we present new results from the Tevatron in this pa- 
per. 

All numbers, unless accompanied by a reference to 
a journal publication, are preliminary. 



II. THE TEVATRON, CDF AND D0 

A. The Tevatron Run II 

The Tevatron in Run II collides protons and an- 
tiprotons at a centre of mass energy of 1.96 TeV 
with a bunch crossing every 396 ns at each interac- 
tion point. Some of the bunches arc by design empty, 
so while the detectors have to be able to cope with 
peak rates of 2.5 MHz, the average collision rate is 
~ 1.7 MHz. Since the start of data taking, the Teva- 
tron has delivered more than 3 fb _1 of integrated lumi- 
nosity at each interaction point, and is now reaching 
peak luminosities of typically 2 • 10 32 cm~ 2 s _1 , with 
the best runs exceeding 2.8 • 10 32 cm^V 1 . Two gen- 
eral purpose detectors take data at the Tevatron, CDF 
and D0. Both collaborations have analysed approxi- 
mately 1 fb _1 of nearly 3fb _1 each has on tape. 

B. The CDF and D0 Detectors 

Both the CDF and the D0 collaboration have a 
strong heavy flavour physics programme, and in the 
last upgrade many features have been added to the de- 
tectors to facilitate this programme. These features 
include precise vertexing, extended \x coverage, and 
sophisticated read-out electronics and triggers. A de- 
scription of the D0 detector can be found at [5j and 
of the CDF detector at Q. 

The same detector characteristics that make D0 
and CDF powerful B physics experiments, also pro- 
vide the basis for a wide-ranging charm physics pro- 
gramme. The feature that makes D0 stand out as 
a heavy flavour experiment particularly is its large \x 
coverage [7|, up to \r\\ < 2. CDF's most distinctive 
detector element for B and charm physics is its fully 
hadronic track trigger Q. 



C. Triggers 

Once every 396 ns there is a bunch crossing at each 
Tevatron interaction point, typically resulting in a 
very busy event; only a few tens of a per mill (ca 
6 ■ 10~ 4 at CDF) of those events can be written to 
tape and analysed in detail, the vast majority will be 
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discarded instantly. The trigger is crucial, and has to 
operate in a challenging environment. There are two 
basic strategies to trigger on heavy flavour events. 

• Trigger on leptons (/x, e). This provides a clean 
signature in an hadronic environment, where 
most tracks are pions. Both experiments use 
this strategy. D0 benefits particularly from its 
large /i coverage, giving large number of semilep- 
tonic and leptonic B and charm events, as well as 
charmonium and bottomium decays, e.g. 180k 
J/ip in the first 0.1 fb" 1 , and 28k B -> J/^K in 
1.6fb _1 (after selection cuts). 

• Trigger on the long lifetimes of B and D mesons. 
This is very challenging in a hadronic environ- 
ment because it requires very fast track and 
impact parameter reconstruction in very busy 
events, in time for the trigger decision. Both 
experiments use this strategy. But only CDF's 
displaced track trigger has enough bandwidth to 
trigger on fully hadronic decays alone, while D0 
have to require an additional lepton in the event 
to keep the trigger rates manageable. This high- 
bandwidth displaced track trigger gives CDF 
unique access to fully hadronic heavy flavour 
decays, such as 13M D° — ► K~ir + events (here, 
and in similar expressions in the rest of the note, 
the charge-conjugate decay is always implied) , 
0.3M D+ -> 0(K+K-)vr+ event as as well as 53k 
B s — ► D^7r T in ~ 1 fb _1 9] (numbers after selec- 
tion cuts). 



III. CROSS SECTIONS 
A. D and B production cross sections 

Amongst the earliest Tevatron Run II results were 
heavy flavour cross section measurements. Figure 
[1] shows CDF's measurement of the prompt differ- 
ential charm production cross section versus pt for 
\rj\ < 1, for D°, D+, D*+ and D s mesons. Using 
only 5.8 pb -1 , less than 0.5% of the presently anal- 
ysed dataset, the statistical error, given by the in- 
ner error bars, is already smaller than the systematic 
uncertainty; the combined statistical and systematic 
uncertainty is given by the outer error bars. The un- 
certainty on the FONLL calculation [ll| is given as 
the grey band, with the central value shown as a solid 
line. While the measured cross sections appear to be 
consistently higher than the FONLL prediction, the 
results are compatible within the (correlated) errors. 

The analyses use the reconstructed impact param- 
eter of the D to distinguish prompt D mesons (with 
zero impact parameter) from those originating from B 
decays (with large impact parameters). 

The same technique of using either impact param- 
eters or decay length has been used to measure inclu- 
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TABLE I: Ratios of measurements and predictions for 
correlated bb cross section measurements, presented by 
F. Happacher at DIS 2006 QJ]. New (2007) results from 
CDF in the m + m~ channel are presented in here. (Table 
taken from [18J], slightly modified.) 



sive B cross sections in B — > D°^ + ^X, B — » D*~ fi + vX. 
and B — > J/ipX decays, where long decay lengths or 
large impact parameters identify D°/i, D*~/i and 3/tp 
originating from B mesons. The inclusive differential 
cross sections measured in the D°/i, D*~/i and the J/tp 
channel are shown in Figures [2] and [3] The results are 
in good agreement with FONLL calculations. The 
preliminary result for the integrated inclusive B cross 
section for 6-hadrons with p t > 9 GeV and \y\ < 0.6 



cr(pp^Hb) = 1.34 //b 

±0.08/xb(stat)+°:ui«b(sys) ± 0.07/ib(BR) 

for the analysis using Hb — > D°(K~7r + )/zX, and 

cr(pp^ H b ) = 1.47 fib 

±0.18 /xb(stat)i°;ig ^b(sys) ± 0.11 ^b(BR) 

for the analysis using H b -> D* + (D°(K-7r+) tt+) /j,X, 
where Hb stands for a generic b hadron. The last 
uncertainty is due to the uncertainty in the branching 
fractions of the specific final states of the D° and D* 
being investigated. The result is in good agreement 
with the FONLL value of 1.39±g;|| |ib [t||. 

CDF performed a measurement of the ex- 
clusive B + — ► J/?/>K + measurement, using 
0.741b" 1 , finding a (p t > 6 GeV, \y\ < 1) 
(2.65±0.12(stai)±0.21(s2/s))^b [i]|. The dif- 
ferential cross sections for the exclusive measurement 
can be seen in Fig [5] together with the B — ► J/ipX 
inclusive results. All measurements disagree signifi- 
cantly with NLO calculations. The agreement with 
FONLL is however very good. 



B. Correlated bb and cc cross sections 

For correlated bb and cc cross sections, i.e. cross 
sections where both the quark and the antiquark are 
within a certain, central rapidity range, higher order 
terms are expected to be smaller and consequently 
NLO calculations are expected to describe the data 
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FIG. 1: Differential charm cross section measured in fully hadronic charm decays using only 5.8 pb -1 at CDF f |lC 
The inner error bars represent the statistical uncertainty, the outer error bars represent the total uncertainty of the 
measurement, and the black line and grey band the central value and uncertainty of the FONLL calculation by [ry. No 
calculation for the D s cross section was available at the time. 
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B -> J/i/>X @| and the FONLL calculation by [T3] - 
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FIG. 2: Comparing inclusive B — > J/^/jX [6| and exclusive 
B + — » JMK + pjjeross section measurements, with Run I 
results Il|], NLO 0, G]| and FONLL [IQJI predictions. 



better [Tj| [ljj]. In the past, Tevatron results on corre- 
lated bb cross sections have been inconclusive. Table U 
shows the status presented by [18] at DIS 2006. Es- 
pecially in the fi + /j," channel, results have been at 
odds with NLO predictions. This year (2007) CDF 
have presented a new measurement of the correlated 
bb and cc cross section in the p + /i~ channel, based on 

pairs with transverse 



0.741b" 1 . The study uses fj, + fi~ 



to bb pairs with px > 2 GeV and a rapidity \y\ < 1.3. 
Prompt fj, are separated from fj, created in charm de- 
cays, and those from jj, created in B decays, using the 
impact parameters of each muon in the /i + /i~ pair. 
The 1-D projection of the impact parameter distribu- 
tion, and fit, is shown in Fig [4] 

The preliminary result for the correlated bb produc- 
tion cross section, where each b quark decays to a p, 
is 



a b^fi,b^fi 



= (1549 ±133) pb. 



and for cc production where each charm quark decays 
to a /i is 



er c _> 



C—*fl,C^*fl 



(624 ± 104) pb 
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FIG. 4: Impact parameter (d) distribution of u, + 1±~ pairs. 
Although only the projection onto 1 dimension is shown, 
the fit is performed on the full 2-D distribution of u, + and 
fi~ impact parameters. The different letter combinations 
that can be found in the legend distinguish muon pairs 
from B ("BB"), from charm ("CC"), and prompt pairs 
("PP") and various mixtures, such as "BP" where one 
muon is prompt and the other from a B decay. 
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FIG. 5: Heavy flavour production mechanisms in pp colli- 
sions. Different mechanisms lead to different kinematics. 



the agreement between measurement and NLO pre- 
diction depends crucially on the choice of e, and no 
firm statements regarding the compatibility between 
CDF's correlated cc cross section measurement and 
NLO calculations can be made until the issue of a 
suitable fragmentation parameter is resolved. 



TABLE II: Measured correlated 66 and cc cross sections in 
0.74 fb" 1 at CDF, divided by NLO predictions, for differ- 
ent assumptions for the Peterson fragmentation parameter 
e M. 



IV. CHARM PAIR CROSS SECTION 



Translating the 6 — ► \x, 6 — > \x cross section to an in- 
clusive 66 cross section, independent of the final state: 

v b - b (PT>6GeV,\y\<l) = 

(1618 ± 148 ± [~ 400 fragmentation]) nb 

where the dominant error comes from the uncertainty 
in the fragmentation function, i.e. the fraction of 6 
quarks that decay to /i. Different values for the Pe- 
terson fragmentation parameter e [20| result in signif- 
icantly different results, especially for the correlated 
cc cross section. This can be seen in Table UH where 
the ratio of the measured cross sections and the NLO 
predictions is given for two values of the Peterson frag- 
mentation parameter e. The value traditionally used, 
e = 0.006j_is obtained from fits to e + ,e _ data [21] . 
However, [16[ point out that these calculations, made 
on the basis of LO parton level cross sections evalu- 
ated with the leading-log approximation of the parton 
shower event generator, cannot be consistently used 
with exact NLO calculations, and find that a more 
suitable value of the Peterson fragmentation param- 
eter would be e = 0.002. For 66 the results for both 
values of e are now in good agreement with NLO pre- 
dictions, in contrast to previous measurements. For cc 



Different production mechanisms for heavy flavour 
quark- antiquark pairs lead to different kinematic dis- 
tributions. The leading production mechanism, and 
their kinematic signature, are depicted in Fig [5] The 
lowest order diagrams ("Flavour Creation") favours 
back-to-back production of the quark-antiquark pair, 
while "Gluon Splitting" favours collinear production. 
Measuring the angular distribution in charm pair pro- 
duction therefore gives clues about the cc production 
mechanism. 

Figure [5] shows the D°, D* _ pair cross section and 
the D + , D*~ pair cross section as a function of the an- 
gle between the two charm mesons, for the kinematic 
range \y D °\ < 1, p%° G [5.5,20]GeV and \y D ' \ < 1, 
p|? * G [5.5, 20] GeV. Using a D* in the reconstruc- 
tion ensures clean data samples. It can be seen that 
in both cases, collinear production is approximately 
of the same size as back-to-back production. The 
figure also shows the result from the Pythia event 
generator, configured to run with leading order ma- 
trix elements plus parton shower (Tune A). The total 
pair production cross sections agree well between the 
Pythia simulation and the data. However, the sim- 
ulation over-estimates back-to-back production and 
under-estimates collinear production. 
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FIG. 6: D , D* pair cross section and D + , D* pair cross 
section as a function of the azimuthal angle between the 
to charm mesons (dots with error bars). The black line 
shows the result from the Pythia event generator [22] , run 
with leadingorder matrix elements + parton shower, tune 
A [H, HI Ell |2<|. The coloured lines show the contri- 
bution to the simulated events from different production 
mechanisms. 



(a) Colour 
Singlet 




(b) QCD with 
higher-order terms 




FIG. 7: "Colour singlet" lowest order QCD (left), and 
QCD with higher order terms (right). The hard scattering 
process is shown in bold. The higher order terms consid- 
ered are those where the "bleaching gluon" is absorbed by 
one of the parton spectators [27] . 



V. CHARMONIUM, BOTTOMIUM 

Producing a colour- neutral J p = 1" state directly 
by gluon-gluon fusion is not possible. The simplest so- 
lution to the problem is to produce a colour-charged 
cc (or bb) pair with gluon fusion and "bleach" it by 
radiating off a gluon. The "colour-singlet" mecha- 



nism, illustrated in Fig[7](a), was first used to describe 
quarkonium production [281 129( , but dramatically fails 
to describe the observed data at the Tevatron, under- 
estimating the 3/ip and ip' production in Run I by 
more than an order of magnitude (see e.g. [3fj). The 
discussion about a solution to this problem is domi- 
nated by two approaches: 

• The "Colour Octet" mechanism proposed by 
[3l| is an effective field theory model based on 
Non-Relativistic QCD (NRQCD). It combines 
results from [32| and elements of the Colour 
Evaporation model (CEM). The CEM was orig- 
inally proposed in [33 . 1341 1 ; m ore recent discus- 
sions can be found in [33, [36J . In the CEM the 
J/ip is essentially formed in a coloured state, so 
no "bleaching gluon" is needed. Then the colour 
"evaporates" in the emission of soft gluons. 

• [22l perform a calculation based on full, rela- 
tivistic QCD, adding higher-order terms to the 
colour singlet term of the type shown in Fig[7]Jb) , 
where the "bleaching gluon" is absorbed by 
a spectator parton. Although each additional 
term is small, there is such a large number of 
such diagrams that the sum of them does indeed 
make a large enough contribution to account for 
the order-of-magnitude difference between the 
observed cross sections and those predicted by 
the colour-singlet model. 

Both NRQCD colour-octet and higher-order perturba- 
tive QCD describe the observed 3/ip and ip' cross sec- 
tions andp t spectra well (see for example [271. [3U |37| ) . 
where the NRQCD approach has a number of ad- 
justable hadronisation parameters that allow a certain 
level of tuning. However, NRQCD makes a firm pre- 
diction that the J/tp should be transversely polarised 
[38j while perturbative QCD predicts a longitudinal 
polarisation of the 3 ftp [27j . 

The ^factorisation approach [H, E3, EH S El 0, 
Ea, Es E3 > can be combined with the colour singlet 
and the colour octet mechanism. In contrast to the 
usual collinear approach, kx factorisation takes the 
non-vanishing transverse momentum of the interact- 
ing gluons into account when calculating the hadronic 
matrix element. Especially when combined with the 
colour-octet mechanism, it describes the production 
cross sections well, although there is a number of ad- 
justable parameters that are not yet fixed by experi- 
mental data [43,Ea|- In contrast to the usual NRQCD 
colour-octet mechanism, kx factorisation predicts a 
longitudinal polarisation of the quarkonium that in- 
creases with pt- 

For a more complete and detailed review of the the- 
ory of quarkonium production at hadron colliders, in- 
cluding comparison to data, see [48| and references 
therein. 
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FIG. 8: Differential cross section of prompt tp(2S) cross 
section in ~ 1 fb -1 , with an earlier result for the J/ip cross 
section |f| superimposed. 



A. The V( 2S ) Cross Section 

The differential cross section of prompt ip(2S) from 
~ 1 fb -1 at CDF Run II and, for comparison the dif- 
ferential cross section of prompt J/ip from 39.7 pb 
at CDF Run II is shown in Fig [8] The prelimi- 
nary result for the integrated cross section of pp — > 
tp(2S)X at y/s = 1.96 TeV, with a subsequent de- 
cay ip(2S) — ► /i + /i~, in the kinematic range |y^(2S)| < 



0.6, p T > 2GeV, is 



a(pp -» ^(2S)X, |^ (2S) | < 0.6, p T > 2 GeV) 
x Br(ip{2S) -> /x+aO 

= (2.60 ± 0.05(stat)t° :ll(syst)) nb 
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FIG. 9: Invariant mass of /i + /j, pairs at CDF (after se- 
lection cuts). CDF finds 0.78M prompt 3/ip in 0.8 fb -1 

sa. 
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FIG. 10: The measured cos 8* distribution in differ- 
ent momentum bins for prompt 3/tp in 0.8 fb - at CDF. 
The pt bins are from left to write, top to bottom: 
[5 - 6], [6 - 7], [7 - 9], [9 - 12], [12 - 17], [17 - 30] GeV, A 
flat distribution indicates no polarisation, a concave (U) 
distribution transverse polarisation and a convex (n) dis- 
tribution longitudinal polarisation. 



B. Measurement of charmonium polarisation 

The J/ip can have three polarisation states, two 
transverse and one longitudinal. In J/ip — > /i + /i~ 
decays, the transverse and longitudinal polarisation 
states can be disentangled by measuring the angle 0* 
between the J/ip and the /i + in the the J/ip restframe. 
It is useful to define the parameter a in terms of the 
cross section of transversely polarised J /ip, ctt, and 
the cross section of longitudinally polarised J/ip, ctl, 
as: 



cr-r — 2a l 
a T + 2a l 



(1) 



If the one longitudinal and the two transverse polari- 
sation states are all equally populated, one would mea- 
sure a = 0; for longitudinal polarisation, a < 0, and 



for transverse polarisation, a > 0. The parameter a 
can be extracted from the distribution of events as a 
function of cos 9* : 



dN 



d{cos0*) 



oc 1 + a cos 



2 a* 



(2) 



In 0.8 fb" 1 , CDF finds 0.8M prompt J/ip -> n + (j,~ 
decays. The mass distribution of the /J, + fi" pairs is 
shown in Fig [5J J/ip originating from B decays are 
rejected using an impact parameter cut. The cos 6* 
distributions are analysed in six bins of different pr, 
as shown in Fig llOl The result in terms of the param- 
eter a as a function of px is shown in Fig[TT](a). The 
plot shows pT-dependent, longitudinal polarisation of 
the prompt 3 ftp. This contradicts NRQCD which pre- 
dicts increasingly transverse polarisation with higher 
momenta; kx factorisation on the other hand appears 
to over-estimate the longitudinal polarisation. The 
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FIG. 11: The polarisation parameter a vs momentum for prompt J/-i/)(a) and prompt i/>(2S) (b) at CDF 49]. Negative 
a correspond to longitudinal polarisation, positive a to transverse polarisation. The NRQCD calculation [3a . [50l . [5l| is 
superimposed in turquoise, the kr factorisation calculation 46] in magenta. 
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FIG. 12: The polarisation parameter a vs momentum for 
prompt J/i/> at CDF Run I and Run II, where the Run II 
data were re-analysed to match the binning in Run I. 




FIG. 13: The fi + fj,~ mass distribution at D0 after selec- 
tion cuts, on the left for all events, on the right only for 
those events selected by the di-muon trigger. 



VI. MEASUREMENT OF THE Y 
POLARISATION 



result also disagrees with Run I results for 110 pb -1 
at a centre-of-mass energy of 1.8 TeV, where evidence 
of positive polarisation of the J ftp was seen. To fa- 
cilitate a comparison between the Run I and Run II 
results, the Run II results have been re-analysed us- 
ing the same binning as in Run I. The comparison is 
shown in Fig[T2l 



The polarisation has also been measured in the the- 
oretically cleaner (no feed down from higher states) tp' 
channel; the result is shown in Fig [11] (b) . The event 
numbers are much lower than for the J/tp but the 
results indicate a trend for the longitudinal polarisa- 
tion fraction to increase with higher px, inconsistent 
with the NRQCD calculation. These results have been 
published in [49J]. 



The theoretical methods describing the production 
and polarisation of T(1S) and T(2S) in p — p collisions 
are equivalent to those for charmonium presented in 
the previous section. 

The D0 experiment measures the polarisation of 
the T(1S) and the T(2S) as a function of the T 
transverse momentum in 1.3 fb -1 of data in Tevatron 
Run II. Data are selected in the di-muon channel. To 
achieve a more reliable estimate of the trigger effi- 
ciency in the cross section calculation, only events se- 
lected by the di-muon trigger are used in the analysis. 
The invariant mass of the (i fi~ pairs near the T mass 
is shown in Fig [13] for all data passing the selection 
cut (a), and those that also were selected by the di- 
muon trigger (b). Figure [l~4l shows the di-muon mass 
spectrum for one bin in | cos9*\ and px, with the fit 
to the data superimposed, and the Gaussians describ- 
ing the individual contribution of the T(1S), T(2S) 
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FIG. 14: The p + At mass distribution at D0 after selec- 
tion cuts, for one bin in cos|0*| and pr, with |cos0*| G 
[0.4,0.5] and p T € [10 GeV, 15GeV]. The plot shows the 
total fit to the data in red for the two fit models briefly 
described in the text, and the Gaussians describing the in- 
dividual contribution of the T(1S), T(2S) and T(3S). The 
background fit is shown as a dotted line. 
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FIG. 16: The polarisation parameter a versus pr mea- 
sured in T(2S) — > n/i events in 1.3 fb -1 of data at D0. 
The data are shown as black circles with error bars. The 
yellow band shows the NRQCD prediction from [53 |. 



[52| is shown as a yellow band, and two limit-cases of 
the ftp-factorisation model [4|| [111 are given as ma- 
genta curves; the flatter of the two magenta curves 
corresponds to A;T-factorisation with the assumption 
of full quark-spin conservation, the steeper curve to 
full quark-spin depolarisation. D0 find clear evidence 
for longitudinal polarisation, in contradiction to the 
NRQCD calculation. The polarisation vs pr in the 
statistically less powerful T(2S) channel, which in- 
cludes contributions from directly produced T(2S) as 
well as feed-down from T(3S) and Xb(2P), is shown 
in Fig [TBJ In contrast to the result for the T(1S), 
the measured T(2S) polarisation is compatible with 
NRQCD. 



FIG. 15: The polarisation parameter a versus pr mea- 
sured in T(1S) — > fifi events in 1.3 fb -1 of data at D0. 
The measured results in eight pr bins are shown as black 
circles with error bars. The yell ow band corresponds the 
NRQCD prediction given in [55q. The magenta lines cor- 
respond to two limit cases of the fcT-factorisation model 
Sali3l- The CDF Run I result [H[ is shown as the green 
triangles with error bars. 



and T(3S). Because of the relatively small number of 
T(3S) events, to ensure a stable fit, the width, rela- 
tive position and the fraction of events in the T(3S) 
peak were taken from MC simulation and fixed in the 
fit. In another fit, the position of the T(3S) peak is 
allowed to float. The difference between the two ap- 
proaches is taken as a systematic error. The T(1S) 
polarisation vs pt is shown in Fig |15l The data repre- 
sent the result for an admixture of directly produced 
T(1S) and T(1S) from other sources, in particular 
from T(2S), T(3S) and Xb- For comparison, the CDF 
Run I result is superimposed in the same plot, as well 
as theory predictions for the polarisation parameter a 
of this admixture. The NRQCD prediction given in 



VII. RELATIVE PRODUCTION CROSS 
SECTION X c2 AND X ci 



The measurement of the Xcj (J = 1,2) cross sec- 
tion is an interesting measurement in its own right, 
as well as important input to J/?/> production mea- 
surements to which it provides an important source 
of feed-down. Experimentally, this measurement has 
always been plagued by the poor mass resolution of 
the reconstructed Xcj which is due to the soft photon 
in the decay chain Xcj - * J/^7- The high luminosity 
at the Tevatron now provides sufficient statistics to 
restrict the analysis only to events where the photon 
undergoes conversions to an e + e~ pair. The momenta 
of the electrons can be measured precisely resulting 
in a far superior energy resolution of the photon than 
would be possible by measuring the photon's energy in 
the calorimeter. This leads to an excellent resolution 
of the reconstructed % c j mass, allowing a separation 
of the J — 0, 1, 2 states. CDF reconstructs Xd in in 
the channel Xci ™ > J/VK/ i / i )7( e+e_ ) within the kine- 
matic range pt(J/'0) € [4, 20] GeV . To select photon 
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FIG. 17: Reconstructed fi + [i (7 — > e + e ) mass, showing 
clearly two distinct mass peaks at the Xci and Xc2 mass. 



conversions, the e + e _ pair is a required to form a well- 
reconstructed vertex a large distance (> 12 cm) from 
the beam, well inside the instrumented region of CDF. 
The spectrum of the (fi + fi~j) invariant mass is shown 
in Fig [TTJ showing two well-separated peaks at the 
Xci mass and the Xci mass. No significant evidence 
for %co production can be seen. After all selection cut 
total number of \d —* 3/ip(jJLfi) 7(e + e~) events found 
in 1.1 fb" 1 is ~7fc. 

The distance between the beamline and the \i + \x~ 
vertex is used to separate the prompt contribution 
from x c j originating from B decays. The measured 
ratio for prompt production is: 



o-{Xc 



o-(Xc 



0.70 ± 0.04(stat) ± 0.03(sys) ± 0.06(BF) 



for p T (J/-0) € [4, 20]GeV. The last contribution to 



the error is due to the uncertainty in the branching 
fraction of Xcj ~~ * J/V'T- This result is at odds with 
the expectation from the colour octet model, which, 
by counting spin-states, predicts a ratio of 5/3. 



VIII. CONCLUSIONS 

The large charm and bottom production cross sec- 
tion at v / s=1.96TeV proton-antiproton collisions, 
combined with the capabilities of the D0 and the CDF 
detector provide the opportunity for new measure- 
ments with unprecedented statistics and precision. 

In this paper we presented Run II results on single 
and correlated open charm and bottom production as 
well as quarkonium production and polarisation. The 
same theoretical framework that managed to describe 
successfully the unexpectedly large charmonium pro- 
duction cross observed at Tevatron Run I now fails to 
account for the significant longitudinal polarisation of 
charmonium and bottomium observed in Run II, nor 
the Xc2 and Xci cross section ratio. Several alternative 
models are being developed, but at the time of writ- 
ing this paper, none has provided a detailed quanti- 
tative post-diction of the charmonium polarisation vs 
Pt that matches the observed data. 

So heavy flavour production is, refreshingly, a field 
that is clearly led by experiment. In response to the 
data from the Tevatron, some of which have been pre- 
sented here, we can look forward to new calculations 
and models, offering new descriptions of the mecha- 
nism of heavy flavour, and in particular quarkonium, 
production. We will also see further new measure- 
ments making use of the large amount of data yet to 
be analysed. Only about 1/3 of the Tevatron data 
taken so far have been used for the results presented 
here, and with the machine going stronger than ever, 
we can expect to at least double the integrated Teva- 
tron luminosity before the end of Run II. 
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